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In their recent Letter Salazar and Toral (ST) study 
numerically a finite Ising chain with non-integrable inter- 
actions 


H/k B T- 


1 ^ t J S^Sj 


2 ^ r 


d+a 


(Si 


±1) 


(1) 


where —d < a < i (like ST, we discuss general di- 
mensionality d). In particular, they explore a presumed 
connection between non-integrable interactions and Tsal- 
lis's non-extensive statistics. We point out that (i) non- 
integrable interactions provide no more motivation for 
Tsallis statistics than do integrable interactions, i.e., 
Gibbs statistics remain meaningful for the non-integrable 
case, and in fact provide a complete and exact treatment; 
and (ii) there are undesirable features of the method ST 
use to regulate the non-integrable interactions. 

ST study a system of finite length L which is simply 
terminated at the ends. Thus the system energy scales as 
E ~ L d+ I CT I (or E ~ L d \aL for a = 0), and an intensive 
energy "density" is obtained from division by a "super- 
extensive volume." We show below that the bulk free- 
energy density obtained from this "non-extensive ther- 
modynamics" depends explicitly on the regulator, so all 
thermodynamics in this boundary-regulated model de- 
pend on boundary effects, and depend on the system 
shape as well for d > 1. A discussion of these issues 
is lacking in Q . 

A preferable model employs periodic boundary condi- 
tions, for which a cutoff in the interaction range must 
be introduced to regulate the energy. Then the bulk 
"non-extensive thermodynamics" depend explicitly on 
the shape of the cutoff, but not on the boundaries or 
the system shape. 

This homogeneous model can be solved exactly, start- 
ing from a modified (nj) with interactions 


J 


R°w{ nj /R) 
(InR^w^/R) 


-d < a < 

(7 = 0, 


(2) 


where the cutoff function w(x) decays at least as fast 
as l/a;' CT ' +E (e > 0) for large x, with w(0) finite. This 
cutoff enables taking the thermodynamic limit. The re- 
maining problem with interactions of range R can be 


mapped exactly to a Kac-potential for cr < by iden- 
tifying 4>{x) — w(x)/x d+a ' , so that the pair interaction 
is R~ d cp{r / R) (we present the case a = elsewhere M). 
One then takes R — > oo, where the lattice becomes negli- 
gible, and recovers the rigorous result for the free-energy 
density (§ 

f/k B T = convex envelope{/°/fcsT - Am 2 } (3) 
where f° is the hard-core free energy, m = L~ d ^ Si, and 


= / TTTfT) io°° w^x^dx -d<a<0 


TT d / 2 w(0)/T(d/2) 


a = 


(4) 


Next, consider reversing the order of limits, with R — > 
oo and finite L. We present only d = 1 for clarity. The 
spin Si interacts with all periodic repeats of Sj, leading 
to a net pair interaction of 


j.. =R a J- w(\kL + nj \/R) 


k=—oo 


\kL 


ll+o- 


(5) 


Remarkably, as R — > oo this sum converges to the con- 
stant value (2/L) w{x)x^~ 1 dx, thus giving pair in- 
teractions that are independent of spatial separation: the 
quintessential mean-field theory. Solving this mean-field 
theory in the thermodynamic limit gives exactly the same 
free energy as before, even for general d, thus demonstrat- 
ing that (Q) is independent of the order of limits. 

Finally, we can treat explicitly the sandwiched case 
L (x R — > oo as well, by use of a hybrid of these two 
methods, obtaining again (|) §. This last limit cor- 
responds directly to "non-extensive thermodynamics," 
since the R a factor in (||) may be interpreted instead as 
the L-dependent temperature employed in pj]. Thus we 
have solved, with standard methods, the homogeneous 
version of "non-extensive thermodynamics," while the 
inhomogeneous version studied in |IJ would result only 
in boundary- dependent modifications to (|^), due to the 
system-shape dependence of the cutoff function. 

In summary, we find non-integrable interactions as 
amenable to Gibbs statistics as integrable interactions, 
leaving the application of alternative methods still with 
burden of motivation. 
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